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ABSTRACT 

Aims. The properties of magnetic fields forming an extended plage region in AR 10953 were investigated. 

Methods. Stokes spectra of the Fe I line pair at 6302 A recorded by the spectropolarimeter aboard the Hinode satellite were inverted 
using the SPINOR code. The code performed a 2D spatially coupled inversion on the Stokes spectra, allowing the retrieval of gradients 
in optical depth within the atmosphere of each pixel, whilst accounting for the effects of the instrument’s PSF. Consequently, no 
magnetic filling factor was needed. 

Results. The inversion results reveal that plage is composed of magnetic flux concentrations (MFCs) with typical field strengths of 
1520 G at log(T) = -0.9 and inclinations of 10° - 15°. The MFCs expand by forming magnetic canopies composed of weaker and 
more inclined magnetic fields. The expansion and average temperature stratification of isolated MFCs can be approximated well with 
an empirical plage thin flux tube model. The highest temperatures of MFCs are located at their edges in all log(T) layers. Whilst the 
plasma inside MFCs is nearly at rest, each is surrounded by a ring of downflows of on average 2.4 km/s at log(T) = 0 and peak 
velocities of up to 10 km/s, which are supersonic. The downflow ring of an MFC weakens and shifts outwards with height, tracing the 
MFC’s expansion. Such downflow rings often harbour magnetic patches of opposite polarity to that of the main MFC with typical field 
strengths below 300 G at log(T) = 0. These opposite polarity patches are situated beneath the canopy of their main MFC. We found 
evidence of a strong broadening of the Stokes profiles in MFCs and particularly in the downflow rings surrounding MFCs (expressed 
by a microturbulence in the inversion). This indicates the presence of strong unresolved velocities. Larger magnetic structures such as 
sunspots cause the field of nearby MFCs to be more inclined. 
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1. Introduction 

In a typical active region on the solar disc three types of features 
can be identified most easily at visible wavelengths: sunspots, 
pores and plage. Whilst sunspots and pores are defined by 
their characteristic darkening of the continuum intensity, plage 
appears brighter than the surrounding quiet Sun mainly in 
spectral lines, or as faculae near the s olar limb in t he continuum. 
It has been known since the work of [Hale ( |1908| ) that sunspots 
and pores harbour magnetic fields with strengths of the order 
of kG. Babcock & Babcock] ( |1955| ) showed that plage, too, is 
associated with magnetic fields, but it was only realised much 
later that it is also predominantly composed of kG magnetic 
features ([Howard & Stenflo| |1972l [Frazier & Stenflo| |1972l 
[Stenflo|1973[); 

The kG magnetic fields, or magnetic flux concentrations 
(MFCs), in plage are often considered to take the form of small 
flux tubes or sheets, and considerable effort has gone into the 
details that determine their structure and dynamics (see the 
review b y Solanki[[1993[). The convective collapse mechani sm 
( [Parker[p~978t [Spruit[[1979| [Grossmann-Doerth et al.|[1998[) is 
thought to concentrate the field on kG values ( [Nagata et al.[ 
[2008t[Danilovic et al.|2010[[Requerey et al.|2014[ ), whereby the 

plasma inside the tube is evacuated and the magnetic field is 
concentrated. This mechanism may not always lead to kG fields, 
however (e.g. [Venkatakrishnan[ [1986} [Solanki e t al.[[1996b[ 
Grossmann-Doerth et al.[[1998[ |Socas-Navarro & Manso Sainzj 


2005[ ). The diameter of an individual kG flux tube is expected 
to be a few 100 km or less, although a lower limit for kG fields 
may exist ( Venkatakrishnan[[ 19 86t [Solanki et al.[[1996a[ ). In the 


internetwork quiet Sun, diameters typically do not exceed 100 
km, necessitating an instrument with an angular resoluti on of 
O'.'15 or better to fully resolve an individual flux tube ([Lagg 
[et al.j2010[ ). 

Owning to the comparatively small lateral size of these flux 
tubes, they ar e commonly trea ted using a thin flux tube model 
( [Spruit[ [T976t [Defouw[[1976[ ), where the lateral variation in 
the atmospheric parameters inside the tube is smaller than the 
pressure scale height. The 3D radiative MHD simulations by 
[Vogler et al.[ ( [2005 [ ) give rise to magnetic concentrations with 
properties that are close to those of the (2nd order) thin-tube 
approximation ( Yelles Chaouche et al.[[200^ . M ore complex 
flux-tube models have also been proposed (see [Zayer et~ar 
( [1989[ ) and references therein). 

Despite their small size, many of the general properties of 
flux tubes residing in plage have nonetheless been determined 
by observations. This has been achieved by analyzing the 
polarisation of the light that is pro duced by the Z eeman effect 
in areas containing magnetic fi eld ([Solanki[[1993[) . Thus, [Rabin 
( [1992[ ); [Zayer et al.[ ( [1989[ ) and [Ruedi et al.[|T992| for example, 
used the deep photospheric infrared Fe I 1.56 pm line to find 
magnetic field strengths of 1400 - 1700 G directly from the 
splitting of this strongly Zeeman sensitive l ine. Fie l d stren gths 
of around 1400 G were a lso obtained by Wiehr[ ( 1978 ) and 
Martinez Fillet et ar| ( [1997[ ) when using lines in the visible, such 


as the 6302 A line pair, whilst values o f only 1000 - 1100 G 


were found by |Stenflo & Harvey[ ( [1985| ) with the 5250 A lines, 
which are formed somewhat higher in the photosphere. Finally, 


the Mg I 12.3 pm lines, used by Zirin & Popp (1989), returned 
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values as low as 200 - 500 G in plage regions, which are fully 
consistent with the kG fields observed in the aforementioned 
lines due to the even greater formation height of the Mg I 12.3 
yum lines, their different response to unresolve d magnetic fields 
and the merging of neighbouring flux tubes dBruls & Solankil 

|T^ . 

The expansion with h eight of MFCs has also been investigated. 
Pietarila et al. ( 2010| ) used SOT/SP images recorded at increas- 


Pillet et al. 


results basec 


( \991) using Milne-Eddington (ME) inversion 


do display internal downflows. 

More recently, observations performed at higher spatial resolu- 
tion u sing the Swedish Solar Telescope (SST; IScharmer et~ar 
2003 1) have further confirmed the basic picture ([Rouppe van der 


Voort et al.|2005| ). |Langangen et al. ( |2007| ) found, by placing a 
slit across a plage-like feature, downflows in the range of 1-3 


ing yu-values and examined the change in the Stokes V signal of 
MFCs in the quiet Sun network. The variations of the Stokes V 
signal across a MFC was, with the help of MHD simulations, 
found to be compatible with a t hin flux tube approximation. 
[Martinez Gonzalez et al.| \20\2) analysed the Stokes V area 
asymmetry across a large net work patch recor ded with IMaX 
([Martinez Pillet et al.j 2011} aboard Sunrise ( [Solanki et al.[ 
[2010t [Barthol et al.[[2Qll[ ) ^ the solar disc centre and found 
that the internal structure of the large network patch was likely 
to be more complex than that of a simple thin flux tube. A 
similar conclusion concerning the internal structure of plages 
was reached by Berger et al. ( 2004| . Rezaei et al. ( 2QQ7[ ) also 


km/s at the edges of the feature. Similarly, Cho et alT] ( [2010[ ) 
observed, using SOT/SP, that pores, too, are surrounded by 
strong downflows in the photosphere. 

The relationship between the magnetic field strengt h and 
continuum intensity of plage was studied extensively by [Kobel 
et al.[ ( [201 1[ ) using ME inversions of SOT/SP data and a clear 


dependence of the continuum intensity on the magnetic field 
strength was found (c.f. Topka et al. 1997[ Lawrence et al. 
[1993 ). Furthermore, the granula r convection in plage areas 
has an abnorma l appearance (e.g . [Title et al.[[1989t [Narayan &| 
Scharmer|2010[ ). [Morinaga et~ar ( [2008[ ) and [Kobel et al.[ ( [2012| ) 


examined the change of the Stokes V area asymmetry of a 
network patch situated at disc centre using SOT/SP and showed 
that it was surrounded by a magnetic canopy. Yelles Chaouche[ 
et al.[ ( [2009[ ) analysed thin flux tubes and sheets produced by 
MHD simulations and concluded that a 2nd order flux tube 
approximation is necessa ry to accurately descr ibe the structure 
of the magnetic features. Solanki et al.[ ( [1999 ) showed that the 
relative expansion of sunspot canopies is close to that of a thin 
flux tube, which could imply that the relative expansion of all 
flux tubes is similar. 

The inclination with respect to the solar surface of MFCs in 
plage was found to be predominantly vertical, with typical 
inclin ations of 10° ( [Topka et al.[[1992[ [Martinez Pillet et ah] 
[1997[ ), which can be attributed to the magnetic buoyancy of 
the flux tubes ( Schussler[[1986[ ), although MFCs with highly 
inclined magnetic fields were also found ([Topka et al.|[1992[ 
[Bemasconi et al.[ [1995t [Martinez Pillet et al.[[1997'| . The 
azimuthal orientation of MFCs was sho wn to have no preferred 
direction ( [Martinez Pillet et al.[ \99l\ and t o form so called 
’azimuth centres’, although Bemasconi et al.[ ( ri995[ ) did find a 
preferred E -W orientation. 

The potential existence of mass motions inside MFCs has been 
fueled by the observation of significant asymmetries in the 
Stokes Q, U, and particularly V p rofiles, in both amplitude 
and area ( [Solanki & Stenflo|[1984[ ). However, [Solanki[ ( [1986[ ) 
showed that within a MFC no stationary mass motions stronger 
than 300 m/s are pre sent. This result was confirmed by [Martini 


concluded that the high spatial density of the kC magnetic fields 
causes a suppression of the convection process. 

As illustrated by the above papers, which are only a small sam¬ 
ple of the rich literature on this topic, there has been significant 
progress in our knowledge of plage in the last two decades. 
Nonetheless, no comprehensive study of plage properties using 
inversions has been published since the work of [Martinez Pillet 
[et al.[ ( p^97] ), which was based on 1" resolution data from the 
Advanced Stokes Polarimeter (ASP). In the following sections 
we aim to both test and expand upon our knowledge of the 
typical characteristics associated with plage using the results 
provided by the recently developed and powerful spatially 
coupled inversion method ( [van Noort|[2012[ ) applied to Hinode 
SOT/SP observations. We concentrate here on the strong-field 
magnetic elements and do not discuss the horizo ntal weak-field 
features also found in active reg ion plage areas ( [Ishikawa et al. 
[20Q8t[Ishikawa & Tsuneta|2QQ9[ ). 


2. Data 

The data set us ed in this investigation w as recorded by the spec 
tropolarimeter ( [Lites & Ichimoto|2013[), which forms part of the 
solar optical telescope (SOT/SP) ([Tsuneta et al.pOOSl'^ 


Suematsu 


on data of the Advanced Stokes Polarimeter (ASP) 
( [Elmore et al.|1992] ). Stokes profiles in plage display a marked 
asymmetry between the areas of their blue and red lobes ( [Stenflo[ 
[et al.|1983| ). This asymmetry is thought to result from the inter¬ 
play between the magn etic element and the convecting plasma 
in which it is immersed ( [Grossmann-Doerth et al. T988[[Solanki[ 
[1989[ ). Following this scenario Briand & Solanki[ ( T998| ) showed 
that low resolution Stokes I and V profiles of the Mg I Z ?2 fine 
can be fitted with a combination of atmospheres representing a 
magnetic flux tube expanding with height, containing no signif¬ 
icant flows, which is surrounded by strong downflows of up to 
5 km/s representing the field-free convecting plasma around it. 
This scenario is generally supported by magneto-hydrodynamic 
(MHD) simulations performed b y e.g. [Deinzer et al.[ ( 1984[); 
Grossmann-Doerth et al.[ ( [1988[ ); Steiner et al.[ ( [1996] ); Vogler[ 


et al. ( 2Q05[ ), although some of the simulated magnetic features 


[et al.[2008] [Ichinioto et al.|2008 Shimizu et al.|2008[ ) aboard the 

Hinode satellite ( Kosugi et al.|2007[ ). The observation was per¬ 
formed on the 30^^ of April 2007, UT 18:35:18 - 19:39:53, using 
the normal observation mode, hence, a total exposure time of 4.8 
s per slit position and an angular resolution of O'/3 was achieved. 
All four Stokes parameters, /, Q, U and V, were recorded at each 
slit position with a noise level of 1 x 10“^ 4. The field of view 
contains a fully developed sunspot of the active region (AR) 
10953 with an extended plage forming region trailing it. During 
the observation the spot was located in the southern hemisphere 
towards the east limb, -190X, -200T, at yu = 0.97 (/i = cos(\6\), 
where 6 is the heliocentric angle). A normalized continuum im¬ 
age of the investigated region used in the inversion is shown in 
Fig.p^ The data were red uced using the standard sp_prep routine 
( [Lites & Ichimoto|2013[ ) from the solar software package. 

3. Inversion Method 

The region of the SOT/SP scan con taining most of the pl age was 
inverted using the SPINOR code ( [Frutiger et al.[[200^ , which 
uses response functions in order to perform a least-squares fit¬ 
ting of the Stokes spectra. It is based upon the STOPRO rou¬ 
tines described by [Solanki et al.[ ( pW7) ). The SPINOR code was 
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extended by van Noort| \20\2) to perform spatially coupled in¬ 
versions using the point-spread-function (PSF) of SOT/SR Such 
spatially coupled inversions have already bee n successfully ap 
plied to Hinode SQT/SP data of sunspots by IRiethmiiller et al. 


(2013 ); |van Noort et al. ( 2Q13| ); [Tiwari et ar] ( 2Q13| ); |Lagg et aT 


( |2014| ). We employ the sa me PSF used by t hese a uthors, which 
is based on the work of |Danilovic et al.| ( |2008| ). The size of 
the inverted area, corresponding to that shown in Fig. is the 
largest that can currently be inverted in a single run by the em¬ 
ployed code due to computer memory limitations. The inver¬ 
sion code allows the recovery of thermal, magnetic and velocity 
gradients with optical depth, among others, which reveal them¬ 
selves by the strengths, shapes and asymmetries present in the 
SQT/SP Stokes profiles (|Solanki|1993[[Stenflo|201Q[|Viticchie &| 


Sanchez Almeida|2011| ). The stratification of each atmospheric 

parameter with optical depth is calculated using a spline interpo¬ 
lation through preset log(T) nodes, where the code can modify a 
pixel’s atmosphere. The resultant full atmosphere is then used to 
solve the radiative transfer equation and the emergent synthetic 
spectra are fitted iteratively by a Levenberg-Marquardt algorithm 
that minimizes the merit function. 

The choice of the log(T) nodes is important for achieving a credi¬ 
ble atmospheric stratification. Three nodes were chosen. A larger 
number of log(T) nodes produced more complex atmospheres at 
the expense of the uniqueness of the solution, while fewer log(T) 
nodes failed to fit the asymmetries in the observed spectra. The 
chosen nodes corresponded to optical depths at log(T) = 0, -0.9 
and -2.3 based on calculated contribution functions of the 630 
nm line pair. The contribution functions were obtained from an 
empirical atmosphere simulating a plage pixel, i.e. containing 
magnetic field of 2000 G at log(T) = 0 and satisfying the thin- 
tube approximation at all heights. We also carried out inversions 
with nodes at slightly different log(T) values, to see if a better 
combination was available, but did not find one for plage. At 
each of the three chosen nodes the temperature, T, magnetic field 
strength, B, inclination relative to the line-of-sight (LOS), y, az¬ 
imuth, line-of-sight velocity, v, and micro turbulence, ^niic^ 
were fitted, leading to 18 free parameters in total. We stress that 
no macro turbulent broadening was allowed and a fixed yu-value 
of yu = 0.97 was assumed during the inversion. The influence 
of straylight from neighbouring pixels is taken into account by 
the PSF and the simultaneous coupled inversion of all pixels. 
Consequently, no magnetic filling factor was introduced in the 
inversion. 

A common problem affecting this inversion process is the possi¬ 
bility that the fitting algorithm finds a solution that corresponds 
to a local minimum. This is particularly so if the initial guess 
atmosphere for a pixel is far from the global minimum. In an 
effort to ensure that the solution for each pixel of the inversion 
corresponds to the global minimum, the inversion process was 
performed a total of four times with each inversion performing 
12 iterations. Save for the initial inversion each successive inver¬ 
sion used the smoothed results of the previous inversion as an 
initial input, thereby ensuring that the initial guess for each pixel 
is closer to the global minimum (under the assumption that the 
inversion does reach the global minimum by itself for the major¬ 
ity of pixels, but runs into danger of falling into a local minimum 
for a minority). After the fourth inversion process the mean;^^ 
value of all the pixels could not be decreased any further, e.g. by 
inverting the scan a successive time. 



Fig. 1. Normalized continuum intensity image of the region in 
which the Stokes profiles were inverted. Pixels where T < 5800 
K are enclosed by the red contour line. The v and y axes indi¬ 
cate the distance to the solar equator and central meridian, re¬ 
spectively. The black box denotes the area taken as a quiet Sun 
reference. 


4. Results 

In this section we describe the various results obtained from the 
inversion. First we give a general overview of the output of the 
inversion followed by subsections dealing with more specific 
points. Figure provides an overview of the continuum inten- 
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Fig. 2. a-c: Magnetic field strength retrieved by the inversion at log(T) = 0, -0.9 and -2.3, from left to right. The colour scale given 
on the right is identical in all the three images, d-f: The line-of-sight inclination of the magnetic field obtained by the inversion 
at log(T) = 0, -0.9 and -2.3, from left to right. All three images have the same colour scale. The black contours in all images 
encompass pixels where T < 5800 K. 


sity, Figs.[2i-c of the magnetic field strength returned by the in¬ 
version, and Figs.[^-f of the inclination, y, of the magnetic field 
vector. All these figures display the entire field of view (FOV) 
to which the inversion code was applied. Fig.l^reveals that part 
of the sunspot’s penumbra as well as pores of various sizes are 
contained in the FOV and had to be excluded from the analysis. 
The sunspot’s penumbra and the largest pores in the image were 


cut out by excluding the lower right hand side of the FOV from 
the analysis. However, many of the pores in the figure are only 
a few pixels in size, illustrated by the contour lines in Figs.[2& 
l^-c and are often entirely embedded within a larger magnetic 
feature. These small pores were removed from the analysis us¬ 
ing a temperature threshold of T < 5800 K at log(T) = 0. Both 
higher and lower temperature thresholds, +150 K, were tested 
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Fig. 4. Vertical slice through a typical MFC. The Y coordinate of this MFC is -151”. a-c\ Magnetic field, LOS inclination and LOS 
velocity from left to right, d-f: Temperature, microturbulence and azimuth from left to right. 


with insignificant effect upon the following results. The 5800 K 
threshold was finally chosen since it corresponds to the lowest 
temperature found in the quiet Sun, (see black box in Fig. [^. 
An (alternative) intensity threshold to remove the pores yielded 
statistically similar results. The area within the black box serves 
as a quiet Sun reference throughout the following sections, since 
it is almost devoid of kG magnetic fields, although hG fields 
are present. The LOS velocities were corrected by subtracting 
an offset of 140 m/s, which was obtained by assuming that the 
pores are at rest on average. 

Figures [^-c reveal that all MFCs expand significantly with 
height, suggesting that many pixels harbour magnetic fields only 
in higher layers of the atmosphere, indicating the presence of 
magnetic canopies. Therefore, all MFC pixels in the inversion 
result were subsequently divided into two populations: core pix¬ 
els and canopy pixels. The core pixels were defined by a pos¬ 
itive magnetic field gradient with optical depth, i.e. a magnetic 
field strength decreasing with height, and an absolute magnetic 
field strength, B, > 1000 G at log(T) = 0. Higher thresholds at 
log(T) = 0 merely reduced the number of selected pixels but did 
not provide results that differ qualitatively from those presented 
here. 

The canopy pixels were defined by a negative magnetic field gra¬ 
dient with optical depth and an absolute magnetic field strength, 
B, above 300 G at log(T) = -2.3. A threshold < 300 G at 
log(T) = -2.3 caused the selection of a large number of pix¬ 
els that were not directly connected to MFCs forming plage re¬ 
gions. These ’extra’ pixels were predominantly associated with 
the sunspot penumbra and canopy, a filament and with weak hor¬ 
izontal magnetic fields found on top of granules in the few quiet 
Sun areas in Fig. This last group is likely related to the weak 
horizontal fields found in plage by |Ishikawa & Tsuneta| ( |2009] ). 
All selected pixels have a Stokes Q,U or V amplitude of at least 
5(T, where cr = 1 x 10“^ A. The location of core and canopy pixels 
using these thresholds is illustrated in Fig. The figure shows 
that the canopy pixels associated with MFCs surround the core 
pixels. The sunspot’s canopy forms an extended ring around it 
and photospheric loops between the sunspot and adjacent, op¬ 
posite polarity pores with field strengths of up to 1000 G at the 
loop apex can be seen. Such a loop structure is located at approx¬ 


imately -215Y, -200T. The sunspot’s canopy extends particu¬ 
larly far, as elongated finger-like structures at -210X, -160F. 
These fingers are presumably very low lying loops connecting 
the spot to MFCs. The clear division between these various mag¬ 
netic structures is only possible by applying the inversion code 
in 2D coupled mode, since the inversion requires no secondary 
atmosphere and/or a filling factor and the remaining single at¬ 
mosphere is given the freedom to differentiate between core and 
canopy fields. 

Most of the MFCs in Figs. |^-f have the same, positive, polarity 
(shown in blue). Only in the lower right hand comer of the field- 
of-view (FOV) can MFCs of negative polarity be found (shown 
in red in Figs.|^-f). The dominant polarity of the MFCs is oppo¬ 
site to that of the sunspot. Between the MFCs of opposite polar¬ 
ities a polarity inversion line (PIL) can be seen, stretching from 
approximately -212Y, -245F to -227X, -210F in Figs.^&[^ 
Ha images, not shown here, indicate the presence of a filament 
along this PIL. The photospheric part of this filament is visible 
in Fig. at log(T) = -2.3 in the form of predominantly horizon¬ 
tal magnetic fields (Fig. of around 350 G. The atmosphere 
below the PIL is almost free of magnetic field (Fig. [^). The 
location of the filament is clearly seen as the elongated canopy¬ 
like structure following the PIL in Fig. For a more detailed 
analysis of th is filament the reader is referred to |Qkamoto et al.| 
( |2008[ |2009| ). Here we can add to their findings that, although 
the filament’s magnetic field reaches down into the photosphere, 
it is largely restricted to layers more than roughly 200 km above 
the solar surface. This geometrical height was obtained from the 
hydrostatic atmospheres returned by SPINOR. The B value of 
350 G in the filament is comparable (within a factor of 2) to 
the field strengths found in AR filaments by Xu et al.| ( |2010| ); 


Kuckein et al. 


( |2012| ); [S^asso et al.] ( |2011| ) in the chromosphere 


sampled by the He I 10830 A triplet. The azimuthal orienta¬ 
tion of the magnetic field within the filament is almost invariant 
across the whole filament. Also, the orientation is not aligned 
with the sunspot’s canopy, but rather almost parallel to the axis 
of the PIL, indicating sheared magnetic fields. This excludes the 
possibility that the field we assign to the filament could merely 
be a low-lying part of the sunspot’s canopy. 

Fig. [^displays a vertical cut through a typical MFC and qual- 


5 








































D. Buehler, et al.: Properties of solar plage 



-250 -240 -230 -220 -210 

X [arcsec] 


Fig. 3. Image showing the location of core and canopy pixels 
using the definition given in Sect. 3. Core pixels are shown in 
orange and canopy pixels are coloured green. The white areas 
contain weak magnetic fields that were not considered to belong 
directly to the MFCs in plage regions. 


itatively illustrates many properties analysed in more detail in 
the following sections. The MFC is composed of nearly vertical 
kG magnetic fields that decrease with height, whilst the MFC 
expands. The apparent asymmetric expansion of the feature at 
-288X arises from the merging of the feature’s canopy with the 
canopy of a nearby MFC. Both the temperature and the micro¬ 
turbulence are enhanced at mid-photospheric layers within the 
MFC, but even more so at the interface between it and the sur¬ 
rounding quiet Sun, where strong downfiows are also present. 
The feature lies between two granules, which can be identified 
in the temperature image at log(T) = 0. The pixel-to-pixel vari- 



Magnetic field [G] 

Fig. 5. Histograms of B values found in MFCs. The three 
coloured histograms are restricted to core pixels, where red 
refers to log(T) = 0, green indicates log(T) = -0.9 and blue 
refers to log(T) = -2.3. The dashed histogram shows the field 
strengths of canopy pixels at log(T) = -2.3. 


ations seen in Fig. are sizeable, but statistically the results are 
quite robust in that they apply to most plage MFCs. 


4.1. Magnetic field strength 

Figures [^-c indicate that the MFCs in plage regions are 
composed of magnetic fields on the order of kG in the lower 
and middle photosphere. This is confirmed by the histograms of 
magnetic field strength in Fig. |5 which are restricted to pixels 
selected using the magnetic field thresholds defined in Sect. 
4. Besides histograms of B of core MFC fields at each optical 
depth, the histogram of the canopy pixels at log(T) = -2.3 
is plotted as well. Histograms of the magnetic field strength 
for the canopy at log(T) = 0 and -0.9 have been omitted as 
at these heights the atmosphere is similar to the quiet Sun or 
contains other, weaker fields that are analysed in Sect. 4.9. 
According to Fig.j^the magnetic field strength at log(T) = -0.9 
has an average value of 1520 G. At this height the two Fe I 
absorption lines show the greatest response to all the fitted 
parameters, making the results from this node the most robust 
and comparable to results obtained from Miln e-Eddington (ME) 
inversions (e.g. Martinez Fillet et al.| \991) of this line pair. 
As expected, the average magnetic field strength in core pixels 
decreases with decreasing optical depth, so whilst at log(T) = 0 
the average field strength is 1660 G, at log(T) = -2.3 the average 
field strength drops to 1180 G. Fig.j^also reveals that the widths 
of the histograms using core pixels decreases with height. At 
log(T) = 0 the FWHM of the histogram is 800 G that then 
subsequently decreases to 580 G at log(T) = -0.9 and to 400 G 
log(T) = -2.3, which is half the value measured at log(T) = 0. 
The comparatively broad distribution at log(T) = 0 appears to 
be intrinsic to the MFCs. Large MFCs display a magnetic field 
gradient across the feature, beginning at one kG at its border 
and rising to over two kG within the space of O'.'5. Smaller 
MFCs, too, often decompose into several smaller features 
when higher magnetic field thresholds are used. However, it 
cannot be completely ruled out, that the field strengths of the 
smallest MFCs are partially underestimated due to the finite 
resolution of SOT/SP. At greater heights neighbouring MFCs 
merge to create a more homogenous magnetic field with a 
smaller lateral gradient in B and appear to loose some of their 
underlying complexity. Nonetheless, differences in the distance 
between neighbouring MFCs still can lead to an inhomogeneous 
magnetic field strength above the merging height of the field 
( Bruls & Solanki||l995| ), which itself strongly depends on this 
distance. 
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BlogH.O [G] 

Fig. 6. Ratio between field strengths at different log(T), J, calcu¬ 
lated using Eq. 3, of core pixels plotted against the local solar 
coordinate corrected LOS magnetic field, B. The dashed line in¬ 
dicates the J of a thin flux tube with .6=2000 G at log(T) = 0. 


The distribution of the canopy pixels indicated by the dashed 
line in Fig. reveals that B in the canopy is generally much 
weaker than in the core pixels. The distribution also has no 
obvious cut-off save for the arbitrary 300 G threshold, suggest¬ 
ing that the MFCs keep expanding with height in directions in 
which they are not hindered by neighbouring magnetic features. 


4.2. Magnetic field gradient 

The change in the peak magnetic field strength in each of the 
coloured histograms in Fig. (see Sect. 4.1) indicates that the 
magnetic field strength of the core pixels decreases with height. 
This was investigated further, first, by correlating the relative de¬ 
crease, d, of the magnetic field for each core pixel using 


6(log(T) = -2.3) 
6(log(T) = 0) ’ 


( 1 ) 


against 6(log(T) = 0), which is displayed in Fig. The dashed 
line in Fig. shows the d predicted by a thin flux-tube model 
with 6=2000 G at log(T) = 0. The thin flux-tube model is identi¬ 
cal to the plage model described by |Solanki & Brigljevic| ( |1992 ). 

Fig. reveals that only magnetic fields above about 2000 G 
at log(T) = 0 have a d close to the plage thin flux-tube model 
( [Solanki & Brigljevic|1992[ ). The greatest limiting factor regard¬ 
ing the calculation of d appears to suffer from a lack of magnetic 
flux conservation between the log(T) = 0 & -2.3 nodes. The flux 
of isolated magnetic features at different logij) was calculated 
by drawing a box around it, which was large enough to easily en¬ 
compass the entire MFC at all heights. Whilst the flux between 
the lower two nodes agreed within 5%, the log(T) = -2.3 node 
consistently contained 20% more flux than the log(T) = 0 node. 
This flux discrepancy between the nodes remained unchanged 
when the flux from a box containing several merged MFCs was 
used. Without this flux discrepancy, a core pixel with 2500 G at 
log(T) = 0 would have a d value of 0.44 instead of 0.55, which 
would bring it close to the predicted value of the plage flux-tube 
model. Between the lower two log(T) nodes, where the flux is 
roughly conserved, the decrease in the magnetic field strength 
with height of such a pixel closely follows the thin flux-tube ap¬ 
proximation. The core pixels with weakest 6 are found at the 
edges of their respective MFC and thus may already be partially 
part of the canopy, due to the limited spatial resolution. This 
would reduce the vertical field strength gradient, thus increasing 
d in particular for core pixels close to one kG. The small oppo¬ 
site polarity patches described in Sect. 4.9 allow the possibility 
of two opposite polarity magnetic fields to exist within a weak 


core pixel at the boundary of the MFC in the lower nodes. The 
Stokes V signals from those two fields would at least partially 
cancel each other, leading to a reduction in the retrieved 6 value 
(and apparent magnetic flux) in the lower two log(T) nodes and 
hence to an increase in d. Another contribution to the mismatch 
in Fig. is that some of the smallest core pixel patches are flux 
tubes which are not fully resolved by Hinode, in particular in 
the lower two layers. The expansion of such an unresolved flux 
tube would then take place primarily within the pixel, leading to 
a nearly unchanging 6 in all three log(T) nodes, i.e. a d close to 
unity. 

The inversion also returns a geometric scale for each pixel. 
However, the inversion process only prescribes hydrostatic equi¬ 
librium within each pixel, but does not impose horizontal pres¬ 
sure balance across pixels. Therefore, each pixel has an indi¬ 
vidual geometric height scale,that can be off-set with respect to 
other pixels. This makes the comparison of gradients in (for ex¬ 
ample) 6 between pixels with very different atmospheres dif¬ 
ficult. However, the core pixels found at the very centre of a 
MFC, with 6 1 2000 G, have very similar atmospheres, allowing 
the estimation of a common gradient in 6. The gradient, A6, in 
the magnetic field of these core pixels is -2.6 + 0.5 G/km be¬ 
tween log(T) = 0 and log(T) = -2.3. As expected from Fig. ^ 
this A6 is smaller than the gradient given by the thin flux tube 
model, which takes a value of A6 = -3.9 G/km over the same 
interval in log(T), but if magnetic flux conservation is imposed 
then the gradient of the inversion agrees with the thin-tube ap¬ 
proximation. 

4.3. Expansion of magnetic features with height 

Fig. 1^ and, in particular. Fig. demonstrate that the MFCs in 
the FOV expand with height. Furthermore, inclination and az¬ 
imuth, reveal that MFCs generally expand in all directions and 
are not subjected to extreme foreshortening effects or deforma¬ 
tions, (see Fig. [^discussed in detail in Sect. 4.6), except due to 
other nearby MFCs (see Sect. 4.7). This raises the question of 
how close this expansion is to that of a model thin flux-tube. 
Several methods were tested to find a robust measure of the 
change in size of a magnetic feature with height. The most obvi¬ 
ous method, the conservation of magnetic flux with height, was 
found to be unreliable to estimate the expansion of the magnetic 
features (see Sect. 4.2). 

The expansion of the MFCs was, therefore, estimated in the fol¬ 
lowing way. At log(T) = 0 the size of a magnetic feature was ar¬ 
bitrarily defined by the number of pixels that harboured a mag¬ 
netic field of at least 900 G. Thresholds above and below this 
value (+200 G) did not significantly alter the results. Then all 
the magnetic field values in the log(T) = 0 image were nor¬ 
malized by the maximum field strength in the feature and the 
ratio, rt, was calculated using Each subsequent 

log(T) layer above log(T) = 0 was in turn normalized by its own 
^maxij) value. The expansion of a magnetic feature could then 
be tracked by the total number of pixels at a given log(T) layer 
where 5 ^^ > rt. This method assumes that the MFCs follow 
a self-similar structure with height. The thin-tube approximation 
as well as some other models (e.g. |Qsherovich et al.|1983| ) fol¬ 
low this principle. Rather than tracking the expansion of a fea¬ 
ture using only the three log(T) nodes returned by the inversion, 
the change of the magnetic field with height was tracked using 
a finer grid of log(T) layers, with a log(T) increment of 0.1. This 
finer log(T) grid was created using the same spline interpolation 
between the three nodes as was used during the fitting by the 
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Relative distance from tube centre R(t)/Ro 



Velocity [km/s] 


Fig. 7. Relative expansion of five isolated magnetic features us¬ 
ing Eq. 2 and shown by the dotted lines. The solid line represents 
the relative expansion of a zeroth order plage and the dashed line 
a network thin fiux tube model. 


inversion procedure (see Fig. |^. The 300 G threshold used to 
select canopy pixels in other parts of this investigation was not 
imposed here in order to avoid an artificial obstruction of the 
expansion. Finally, the relative expansion of a feature was calcu¬ 
lated using 


Rir) 

Ro 


A(t) 
Ao ’ 


( 2 ) 


where R and A are the radius and area, respectively, of the fiux 
tube at optical depth r, Rq and Aq at log(T) = 0. 

Five isolated magnetic features were selected from within the 
field of view. The number of selected features is small since most 
magnetic features have merged with other features at log(T) = 
-2.3, as demonstrated in Figs. [^&[^ The R(t)/Ro of the five 
selected features are represented in Fig. [7] by the dotted curves, 
while the solid line shows the relative radius of the 0^^ order 
thin fiux-tube plage model of |Solanki & Brigljevic| ( |1992| ). All 
the dotted lines in Fig. |7] follow the expansion predicted by the 
plage model reasonably we ll. Interestingly , the thin fiux-tube 
model for the solar network ( |Solanki|1986| ), dashed line in Fig. 
|7] did not fit the expansion of the observed MFCs as well as the 
plage model {solid line in Fig. [^. The reduced relative expan¬ 
sion of the selected features above log(T) = -2 when compared 
to the model may be an indication of the merging of features 
limiting the expansion at those heights. Another possibility is 
that a zeroth order model is not sufficient to describe the expan- 
sion of the selected fe atures, especially in higher layers ( |Yelles| 
Chaouche et al.|20d9| ), since the lateral variation of the magnetic 


field within the tube is no longer negligible in higher order fiux 
tube models ( [Pneuman et al.||1986| ). Given that the majority of 
MFCs merge with nearby MFCs, it follows that the majority of 
MFCs are expected to depart from the expansion displayed by 
isolated MFCs. 


4.4. Velocities 

Figure [8 i-c displays the LOS velocities retrieved by the inver¬ 
sion at the three log(T) nodes. The small FOV for this figure, 
representing a typical plage region, was chosen to better high¬ 
light the striking differences between velocities in field-free and 
kC regions and the unusual velocities recorded at the interface 
between these two regions. The core pixels in the images are 
enclosed by the thin black contour lines. Outside the areas har¬ 
bouring core pixels, the typical granular convection pattern can 
be seen at log(T) = 0 & -0.9. The LOS velocities in the top 
node outline only traces of the stronger granules and display 
some similarities with chromospheric observations, albeit with 


Fig. 9. Histograms of the LOS velocities found in MFCs. The 
coloured histograms were obtained for core pixels, where red 
refers to log(T) = 0, green to log(T) = -0.9 and blue refers to 
log(T) = -2.3. The dashed histogram shows the los velocities of 
canopy pixels at log(T) = -2.3. 


smaller velocity amplitudes. The dotted contour lines outline the 
canopy fields. 

The plasma within the majority of core pixels is nearly at rest in 
all log(T) nodes, as Figs.[^-c qualitatively indicate. Histograms 
of the LOS velocities of these pixels displayed in Fig. [^support 
this assertion. The mean and median velocities of these core pix¬ 
els are 0.8 km/s and 0.6 km/s, respectively, at log(T) = 0. They 
drop to 0.2 km/s and 0.2 km/s at log(T) = -0.9, and to 0.1 km/s 
and 0.0 km/s at log(T) = -2.3. The red histogram in Fig.[^ cor¬ 
responding to the log(T) = 0 layer, also contains a significant 
fraction of pixels with downflows larger than 1 km/s. The tail of 
faster downflows is mainly responsible for the larger than aver¬ 
age velocity at log(T) = 0. The individual MFCs, one of which 
is displayed in Fig. were inspected further to determine the 
location and nature of these fast downflows at logr = 0. Fig. 
[TO^ reveals the core pixels of a MFC, enclosed by the black con¬ 
tour line, to be surrounded by a ring of strong downflows. Other 
MFCs show similar downflow rings at log(T) = 0 and can often 
be identified based on such a ring alone. Downflows that exceed 
1 km/s are never found within a MFC, but occasionally a core 
pixel located at the edge of a MFC can coincide with a strong 
downflow, giving rise to the tail of strong downflows seen in the 
red histogram in Fig. Furthermore, the velocities within the 
rings are not of uniform magnitude. Often some portions of a 
ring show very fast downflows, up to 10 km/s, whilst others can 
harbour flows of barely 1 km/s. The portions featuring the fastest 
downflows do not have a positional preference with respect to its 
MFC, of which Fig.[^is one example, which precludes that the 
magnitude of the downflow is affected by the viewing geometry. 
The downflow ring seen at log(T) = 0 in Fig. [T^ is still visible 
at log(T) = -0.9, in Fig.p^. However, the pixels with the fastest 
downflows in the ring seem to be located further away from the 
core pixels in this layer, when compared to Fig. 10 1 . It appears 
that the downflow ring shifts outwards as the magnetic field of 
the MFC expands with height (see Sect. 4.3). Also, the ring ap¬ 
pears to be wider at this height. At log(T) = -2.3 the ring can no 
longer be identified. 

A quantitative picture of the LOS velocities within these rings 
was gained by analyzing the pixels, which directly adjoin the 
core pixels. Only the lower two layers were analysed and are 
displayed in Fig. since the rings are no longer present at 
log(T) = -2.3. LOS velocities of up to 10 km/s were found 
within these rings at log(T) = 0 and the mean and median values 
for the corresponding histogram are 2.44 km/s and 2.16 km/s, 
respectively. For comparison, the fastest downflow in the quiet 
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Fig. 8. a- c: Line-of-sight (LOS) velocities retrieved by the inversion at log(T) = 0, -0.9 and -2.3, from left to right. The thin black 
contours outline core pixels and the black areas at log(T) = 0 mark supersonic velocities. The dotted lines (in panel c) display 
canopy pixels. 
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Fig. 10. Same as Fig.[^ but for a single MFC. 
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Fig. 12. Same as Fig.[^ but for a small part of the full FOV, chosen to reveal a downflow plume around a MFC. The solid contour 
line bounds core pixels and the dotted line (in panel c) canopy pixels. The arrows point to the location of the plume at each height. 


Sun region, see the black box in Fig. is only 6 km/s. At 
log(T) = -0.9 the rings no longer contain downflow velocities 
faster than those found in the quiet Sun at the same log r layer, 
but the histogram of the velocities in the ring still has a mean 
LOS velocity of 0.84 km/s and a median of 0.77 km/s. We further 
investigated whether the fast downflow velocities at log(T) = 0 
in these rings attain supersonic values. Since the SPINOR code 
calculates a full atmosphere, including density and pressure, for 
each pixel, we were able to directly calculate a local sound speed 
for each pixel over it s entire log(T) range, using the same ap¬ 
proach as |Lagg et al. ( |2014| ). The adiabatic index, also needed 
to calculate the sound speed, was acquired from look-up tables 
produced by the MURaM MHD simulation code ( [Vogler et al.| 


|2005| ). Supersonic velocities were found in pixels with down¬ 
flows exceeding 8 km/s at log(T) = 0 and the fastest downflows, 
reaching 10 km/s, have a Mach number of 1.25. Higher log(T) 
layers did not show any supersonic velocities in any of the pix¬ 
els. Since the highest speed found in the quiet Sun reach up to 
6 km/s, no supersonic velocities were consequently found in the 
quiet Sun. Furthermore, we determined that 2.5% of a MFC’s 
downflow ring contain supersonic velocities. Pixels containing 
supersonic velocities are coloured black in Fig.[^ 

Whilst the downflow rings seen at log(T) = 0 & - 0.9 are gen¬ 
erally not traceable at log(T) = -2.3, many MFCs additionally 
have downflows in the form of a plume-like feature, which can 
be traced through all three logij) layers. An example of such a 


9 







































































D. Buehler, et al.: Properties of solar plage 



Velocity [km/s] 

Fig. 11. Histograms of the LOS velocities of pixels immediately 
surrounding core pixels, where red refers to log(T) = 0, green to 
log(T) = -0.9. The dotted histograms display LOS velocities in 
the quiet Sun. 


plume-like feature is displayed in Fig. T he plume lies just 
outside the core pixels at log(T) = 0 (Fig. 1 12^), but is consider¬ 
ably further away, at the boundary of the canopy at log(T) = -2.3 
(Fig.p^), and appears to trace the expansion of the MFC. It also 
increases in size with height. At all heights the LOS velocities 
of the feature are high when compared to their immediate sur¬ 
roundings, but only at log(T) = 0 are the velocities in the feature 
higher than can be found in the quiet Sun. As with the downflow 
rings the velocities progressively increase with depth. 


4.5. Temperature 


The temperatures at each of the three log(T) nodes is displayed 
in Fig. for the same FOV as in Fig.[^ Fig.[T^ corresponds 
to the temperature at log(T) = 0 and exhibits the familiar 
granulation pattern. The positions of core pixels are revealed by 
the black contour line in the images and show that they (the core 
pixels) are predominantly found within the comparatively cool 
intergranular lanes. Figures [T^ &[T^ display the temperature 
at log(T) = -0.9 and -2.3, respectively. Both images indicate 
the high temperatures within MFCs when compared to the 
quiet Sun at these heights. Furthermore, the temperature map at 
log(T) = -2.3 not only reveals the comparatively hot MFCs but 
also a reversed granulation pattern. 

The higher temperatures within MFCs at log(T) = -0.9 and 
-2.3, compared to the quiet Sun, are illustrated more quan¬ 
titatively by the histograms in Fig. lA At both those layers 
the average temperature is around 300 K higher within core 
pixels, with average temperatures of 5690 K and 5070 K at 
log(T) = -0.9 and -2.3, respectively, than in quiet Sun pixels, 
where the average temperatures at the same log(T) heights are 
5290 K and 4780 K, respectively. Fig.p^ further demonstrates 
that the average temperature in canopy pixels at log(T) = -2.3 
is only slightly lower than the temperatures of core pixels at the 
same height, with a mean of 5000 K. Only at log(T) = 0 is the 
average quiet Sun temperature higher, at 6410 K, than in the 
core pixels, which have a mean temperature of 6270 K. Quiet 
Sun pixels harbouring downflows {dotted red histogram in Fig. 
[T4^ ) have a slightly lower mean temperature at 6240 K than 
MFCs, which are also located predominantly in downflowing 
regions. The MFC temperature histograms at log(T) = 0 have 
been artiflcially curtailed by the 5800 K threshold imposed at 
the beginning of the investigation to remove pores. However, 
the bulk of plage pixels is well above this threshold. Quiet Sun 
pixels in upflowing regions at log(T) = 0 {dotted blue histogram 


in Fig. Mi) have a mean temperature of 6590 K, well above the 
values of MFCs. 

The temperature gradient within core pixels was studied further, 
first, by taking the ratio of the log(T) = 0 and -2.3 temperatures. 
A histogram of these ratios is seen in the left panel of Fig. 
P?| The average ratio for core pixels is 0.81 + 0.02, which 
demonstrates that the majority of core pixels have a very similar 
temperature stratification. The ratios were then compared to the 
te mperature ratio obtained from a plage flux tube model derived 
by Solanki & Brigljevic ( 1992|), which is shown by the dotted 
line in the left panel in Fig. 15 The temperature ratio of the 
model, which has a ratio of 0.79, agrees reasonably well with 
the inversion results. The thin flux-tube network model ( [Solanki | 
[1986 ), has a ratio of 0.7 between the same log(T) heights and 
lies outside the histogram. The temperature stratification of 
MFCs studied in this investigation significantly deviate from the 
network model’s prediction, which is not surprising given that 
we are studying strong plage. 

The temperature stratification of various models is depicted in 
the right panel of Fig. The dotted and dot-dashed curves 
represent the plage and network flux-tube models ( [Solanki 


1986[ Solanki & Brigljevic [1992[ ), while the solid line in the 


same panel shows the typical temperature stratification obtained 
from core pixels by the inversion. Whilst quantitatively the 
model and inversion results agree quite well, in particular at the 
three log(T) nodes, qualitatively there is an important difference. 
The temperature stratification of the model has a notable bend 
between log(T) = -1 and -1.5, which is entirely absent from the 
inversion result. However, such a bend can only be reproduced 
by the inversion by employing at least four nodes, which would 
in turn introduce too many free parameters, when inverting only 
2 spectral lines. 

Whilst on average the temperature stratification of a MFC 
follows that a thin flux tube model, a closer inspection of the 
MFC’s cross-sections reveals that the temperature within a MFC 
is not uniformly distributed at all. The highest temperatures 
appear to be preferentially located at the edges of the MFCs. 
At log(T) = 0, Fig. [^, the temperature gradients across a 
MFC are strongest, whilst higher layers display an ever more 
uniform distribution of temperatures. Nevertheless, even at 
log(T) = -2.3, Fig. [^, some areas within the MFCs have 
an enhanced temperature when compared to their immediate 
surroundings. These enhanced temperature areas can usually be 
traced through all three layers and become smaller in size in 
deeper layers, e.g. at -242.5X and -235T or at -238.5X and 
-240F. The white areas in Fig. [T3^ indicate the pixels contain¬ 
ing supersonic velocities (see Sect. 4.4). Many examples can be 
found in Fig. of a match between the location of supersonic 
velocities and a nearby (1-2 pixels away, but always within the 
MFC) local temperature enhancement. However, there is no 
one-to-one relationship between the two. At numerous locations 
throughout the plage region, e.g. at -238X and -237.5F in 
Fig. there is a clear local enhanced of the temperature 
across all three log{T) layers, but no nearby supersonic velocity. 
Furthermore, no linear relationship seems to exist between the 
magnitude of the temperature enhancement (in any layer) and 
the magnitude of the nearby supersonic downflow. 


4.6. Inclination & Azimuth 

The inclinations of the magnetic held plotted in Figs. are 
inclinations in the observer’s frame of reference. Due to the 
small heliocentric angle {{6) = 13^) a qualitative picture of the 
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Fig. 13. a-c: Images of the temperature at log(T) = 0, -0.9 and -2.3, respectively. The black contour lines encompass core pixels 
and the white areas at logir) = 0 denote pixels containing also supersonic velocities. The dotted lines (in panel c) display canopy 
pixels. 
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Fig. 14. a\ Temperature histograms at log(T) = 0 of core plage pixels, solid line, and the quiet Sun {dotted lines). The quiet 
Sun temperatures have been further divided into temperature histograms corresponding to downflowing, dotted red, and upflowing, 
dotted blue, regions, b: Temperature histograms at log(T) = -0.9 of core plage pixels, solid, and the quiet Sun, dotted, c: Temperature 
histograms at log(T) = -2.3 of core plage pixels, solid, and the quiet Sun, dotted. The dashed histogram displays the temperatures 
in canopy pixels. 


inclinations of plage magnetic features can still be gained from 
those figures. However, a conversion of these inclinations to 
local solar coordinates is necessary to find the inclination of the 
magnetic fields with respect to the solar surface. 

The conversion of the inclinations and azimuths retrieved by 
the inversion to local solar coordinates is not straightforward. 
Whilst the inclination of the magnetic field is uniquely defined 
the azimuth has an inherent 180° ambiguity ( [Unno| p^56| ). 
Therefore, when converting to local solar coordinates one is 
forced to choose between one of two possible solutions for the 
magnetic field vector. Many codes, requiring various amounts of 
manual input, have been developed to so l ve the 180° ambiguity. 
The re ader is directed to [Metcalf et al.| ( |2006| ) and [Leka et aT] 
( |2009j ) for overviews. An additional challenge facing these 
codes is that they generally use the output of an ME inversion as 


an input. The output of the ME inversion does not contain any 
information on the change of the magnetic field vector over the 
formation height of the inverted absorption line. The SPINOR 
code, however, provides this information, which in turn allows 
the canopies of magnetic features to be identified as is shown 
in Eig.|^ The azimuths retrieved by the inversion were spatially 
very smooth in all the nodes, indicating that the azimuth was 
well defined in the regions containing magnetic fields. 

The canopy pixels shown in Eig.[^form continuous rings around 
the various cores. By assuming that the magnetic field in each 
canopy pixel originates from the largest patch of core pixels in 
its immediate vicinity, the direction of the magnetic field vector 
of each canopy pixel can be determined unambiguously as long 
as the polarity of the corresponding core patch is known. The 
polarity of a patch of core pixels can be determined easily from 
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Ratio 



Fig. 15. Left: Distribution of temperature ratios of the log(T) = 0 and -2.3 layers of all core pixels. The dotted line shows the 
same ratio obtained from the plage flux tube model of |Solanki & Brigljevic| ( |1992] ). Right. The solid line represents the temperature 
stratification of a typical core pixel obtained from the inversion. The dotted line follows the temperature stratification of an ideal 
plage fiux tube model and the dot-dashed line indicates the temperature stratification of the network fiux tube model. The dashed 
line depicts the temperature stratification of the HSRA. 
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Fig. 16. a-c\ Ambiguity resolved azimuths in local solar coordinates, O, at log(T) = 0, -0.9 and -2.3 from left to right. North is up 
and corresponds to an angle of 90°. West is to the right and corresponds to an angle of 0°. d-f\ The local solar coordinate corrected 
inclinations, T, of the magnetic field after the azimuth ambiguity resolution at log(T) = 0, -0.9 and -2.3 from left to right. The 
inclinations and azimuths of pixels with 5j300 G at log(T) = -2.3 and 5j700 G at log(T) = 0 & -0.9 are shown in white. 


their Stokes V spectra. In essence, the magnetic field vector of 
a canopy pixel points towards a patch of core pixels if it has 
a negative polarity and, in turn, away from it if the patch of 
core pixels has a positive polarity. With the help of this process 
we were able to resolve the azimuth ambiguity of the canopy 
pixels without having to make any further assumptions on the 
properties of the magnetic field. This process was repeated for 
each canopy pixel individually. Since the canopies of the various 
plage features are greatest at log(T) = -2.3, the magnetic field 
vector of the canopy pixels was determined at this log(T) height. 
Once the magnetic field vectors of the canopy pixels were 
determined, the vectors of the core could be obtained using 


an acute-angle method. This method works by performing dot 
products, using the two possible vectors in an undetermined 
pixel, with those surrounding pixels whose vector was already 
determined. The dot products corresponding to each of the two 
possible solutions were then summed. The vector associated 
with the smallest sum was subsequently selected as the correct 
vector for that pixel. The pixels that were surrounded by the 
largest number of already determined field vectors were given 
preference. 

Once the magnetic field vectors of both the canopy and core 
pixels at log(T) = -2.3 were known, the vectors at log(T) = 0 
& -0.9 could also be determined. The now known vector at 
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Fig. 17. Histograms of magnetic field inclination of MFCs rel¬ 
ative to the solar surface normal, F. The three coloured his¬ 
tograms were obtained using core pixels, where red refers to 
log(T) = 0, green shows to log(T) = -0.9 and blue refers to 
log(T) = -2.3. The dashed histogram depicts F of canopy pixels 
at log(T) = -2.3. 



Azimuth, O [degrees] 

Fig. 18. Histograms of O found in plages. The three coloured 
histograms were obtained using core pixels, where red refers to 
log(T) = 0, green shows to log(T) = -0.9 and blue refers to 
log(T) = -2.3. The dashed histogram shows O of canopy pix¬ 
els at log(T) = -2.3. The dotted line represents a homogeneous 
distribution. 


log(T) = -2.3 in each pixels was used to perform dot products 
with the two possible vectors in the next lower log(T) layer. The 
vector with the smaller dot product was subsequently chosen 
as the correct magnetic field vector. The 180° ambiguity of the 
magnetic field vector at log(T) = 0 & -0.9 was removed for 
only those pixels where 5^700 G in either layer. The resolution 
of all the vectors is entirely automatic and only the definition 
of the core and canopy pixels for the initial input is manual, 
but followed the definition given in Sect. 3. Also, no smoothing 
of the azimuths is performed at any point. The converted 
inclinations and azimuths in local solar coordinates after the 


resolution of the 180° ambiguity are plotted in Fig.p^ 

Figures [T^-c show the resolved azimuths, O, at all three nodes. 
Several ’azimuth centres’ ( [Martinez Fillet et~ST]|1997| ) can be 
readily identified, in particular at log(T) = -2.3. In combination 
with Figures [T^-f it can be seen that most of the ’azimuth 
centres’ have vertical fields in their cores that become more 
horizontal towards the edges of a MFC. ’Azimuth centres’ tend 
to be either relatively isolated features or large ones. Most of 
the MFCs tend to be elongated with the field expanding roughly 
perpendicular away from the long axis of the structure over 
most of its length and directed radially away at the ends. 

A more quantitive picture of the general inclinations of the 
core pixels can be obtained through their histograms, depicted 
in Fig. [T^ The distributions of the inclinations have their peak 
between 10° and 15° for all log(T) nodes. The mean inclination 
for each log(T) layer is 22°, 18° and 23° with decreasing optical 
depth. The median values have a similar progression with 
optical depth, taking values of 19°, 17° and 21° respectively. 
Fig. also reveals that the canopy pixels are significantly more 
horizontal, with a peak at 25° and a very extended tail reaching 
up to 90°. The mean inclination for the canopy fields is 39° 
with a median of 36°, which demonstrates quantitatively the 
more horizontal nature of the canopy when compared to the 
core fields. The largest inclinations are found at the edges of the 
canopies as expected for an expanding fiux tube or fiux sheet. 
The histograms of O of core pixels are depicted in Fig.[^ None 
of the four distributions are homogeneous and show a consistent 
under-representation of the W and partly the N directions. These 
two directions are, however, expected to be under-represented 
due to the viewing geometry, as the region was located in the 
SE dii the time of the observation. The azimuth distributions 
from MFCs found at the northern edge of the field of view show 
a more homogeneous distribution, as expected. 

Fig. [T^ also shows that the peak of each azimuth distribution 


is shifted with respected to other distributions, suggesting that 
the direction of the magnetic field vector of individual pixels in 
MFCs appears to rotate with height. Several tests were carried 
to determine the nature of this rotation, after which a solar 
origin as well as an inversion based error seem unlikely. Several 
instrumental effects such as cross-talks or differences in the 
spectral dispersion between the Stokes parameters were found 
to be capable in causing the observed rotation. However, a more 
in depth investigation regarding this matter is required. 


4.7. Effect of the sunspot 


The majority of the MFCs show no obvious and conclusive 
influence of the nearby sunspot’s magnetic field, which stretches 
well beyond the spot’s visible boundary in the form of a 
low-lying magnetic canopy, as can be seen in Fig. Similar 
extended sunspot canopies were already found in earlier studies 
(e.g. |Giovanelli|198^|Solanki et al.|19^|1994| ). A few MFCs 
in the FOV, however, are noticeably influenced by the sunspot. 
The most striking of these features is located at -217X, -125F 
in Figs.|^&[^ where an extensive loop system can be identified. 
These loops have horizontal fields and connect several pores 
with positive polarity to the negative polarity sunspot. The 
magnetic field strengths found in this loop system can reach 
values as high as 1000 G at log(T) = -2.3 in a few places 
and Fig. indicates that they are suspended above relatively 
field-free gas, since they are identified as canopy pixels. 

MFCs close to the sunspot also posses highly deformed 
canopies, which are elongated towards the spot if the MFC has 
the opposite polarity of the spot. Such MFCs can be seen at 
-220X,-170Y and histograms of their inclinations and azimuths 
are presented in Fig. [T^ The field in the core pixels of these 
MFCs is more inclined than on average; compare with Fig. 
[TtI The mean inclinations of the field at the three layers from 
log(T) = 0 to log(T) = -2.3 are 31°, 32° and 35°, respectively. 
These average inclinations are about 10° larger than for MFCs 
found further away from the spot. In particular the inclinations 
of the canopy pixels in Fig. [T^ reveal the effect of the sunspot 
upon these magnetic fields even more strongly. The mean 
inclination of the canopy fields is 56° and the median value is 
57°, which is more than 15° larger than the average in Fig.[T7| 
The azimuth distributions in Fig. [^clearly display the influence 
of the sunspot as all the distributions both from the core and 
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Inclination, r [degrees] 



Fig. 19. Left: Histograms of F found in plage around -220X, -170F. The three coloured histograms were obtained using core pixels, 
where red refers to log(T) = 0, green shows to log(T) = -0.9 and blue refers to log(T) = -2.3. The dashed histogram represents 6 of 
canopy pixels at log(T) = -2.3. Right: Histograms of O found in plages around -220X, -170F. The three colours and the dashed 
line have the same significance as in the graph on the left in this figure. The dotted line represents a homogeneous distribution. 


canopy pixels show a clear preferred orientation towards the 
spot. It is likely that the canopies of these MFCs interact with 
the sunspot’s canopy ( [Solanki et al.|1996^ , even if the fields of 
both magnetic structures do not appear to be directly connected 
over the log(T) range considered in this investigation. 

From the examples given in this section the following picture 
emerges. Those magnetic features that are in the immediate 
vicinity of the sunspot, i.e. if there are no further kG magnetic 
fields in between them and the spot, show a clear deformation of 
their canopy and have inclination and azimuth distributions that 
either predominantly point towards or away from the spot, de¬ 
pending on the polarity. The spot’s infiuence on the orientation 
of magnetic features could be observed up to 20” away from 
the sunspot’s outer penumbral boundary, provided there were 
no other magnetic features in between. Any magnetic feature 
situated behind this ’first row’ of kG features is effectively 
shielded from the spot in the photosphere and then behaves 
closer to an isolated magnetic feature, affected only by its 
nearest neighbours. 


4.8. Weak opposite polarity fields next to MFCs 


Figures [2Q^-c show the LOS inclinations of the magnetic field 
over a small FOV (same as Figs.[^[^&|^. In this blow-up the 
fields in the MFCs {black contour lines) display mainly vertical 
orientations, shown in blue, but a closer inspection of the images 
corresponding to the log(T) = 0 & -0.9 layers exhibit that many 
MFCs are adjoint by magnetic patches with an inclination op¬ 
posite to the MFC, as can be seen from their red colour, e.g. at 
-237.5X and -238.5F. Further examination of Fig, [^reveals 
that these small opposite polarity patches are hidden beneath the 
canopy of the nearest MFC. As a result the small red patches are 
absent in Fig. [20]:. The vast majority of these small opposite po¬ 
larity patches have B\ 100 G at log(T) = 0 & - 0.9, whereas the 
canopies above them have 5^300 G. This means that the Stokes 
spectra of these pixels are dominated by the canopy fields. Only 
the deconvolved Stokes profiles returned by the inversion, as dis¬ 
played in Fig.[^ show a small polarity reversal in the wings of 
the Stokes V profile. 

Nonetheless a number of pixels in these small opposite polarity 
patches have field strengths in excess of 100 G at log(T) = 0. 
These pixels were analysed further by producing histograms of 
their field strengths and LOS velocities, which are depicted in 
Fig. [22^ & b, respectively. At log(T) = 0 the mean and me¬ 
dian of the histogram are only 280 G and 190 G, respectively. 



6301.0 6301.5 6302.0 6302.5 6303.0 

deconvolved Wavelength [A] 

- fitted 

- observed 

Fig. 21. Stokes spectra of a typical canopy pixel harbouring a 
small opposite polarity magnetic field at log(T) = 0 and -0.9. 
The black spectra correspond to the original SOT/SP obser¬ 
vation, the red curves display the spectra fitted by the inver¬ 
sion, and the resultant spectra after the spatial deconvolution are 
coloured green. 


These values drop to 160 G and 70 G at log(T) = -0.9 before 
rising to the considerably higher strengths of the canopy fields 
at log(T) = -2.3. The histograms in Fig. [2^, corresponding to 
the lower two log(T) layers, demonstrate that these weak oppo¬ 
site polarity fields are predominantly located in the strong down¬ 
flows surrounding the MFCs. The mean and median velocities of 
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-242 -240 -238 -236 -242 -240 -238 -236 -242 -240 -238 -236 

X [arcsec] x [arcsec] x [arcsec] 

Fig. 20. a-c\ LOS inclination, y, at log(T) = 0, -0.9 and -2.3, respectively. The black contour lines encompass core pixels. The 
arrows point to the location of the weak opposite polarity at each height. 




Magnetic field [G] Velocity [km/s] 

Fig. 22. a: Histograms of the magnetic field strengths in pixels featuring weak opposite polarities. The colours red, green, blue refer 
to log(T) = 0, -0.9 & - 2.3, respectively. b\ Histograms of the LOS velocities in the same pixels. The colour code is identical. 


these pixels at log(T) = 0 are 1.6 km/s and 1.3 km/s, respectively. 
At log(T) = -0.9 the mean and median velocities are 1.3 km/s 
and 1.2 km/s. The velocities at log(T) = -2.3 are, as expected, 
the same on average to those seen in the dashed histogram in 
Fig. 1^ since at this height the LOS velocity inside the MFC is 
sampled and, hence, mostly weak fiows are seen. The fiux stored 
in these opposite polarity patches is tyically below 1 % of the fiux 
of their parent MFCs. 


4.9. Microturbulence 


The inclusion of micro- and macro-turbulent velocities is com¬ 
mon when fitting photospheric absor ption lines and is an indica¬ 
tion of unresolved fin e structure (e.g. |Lites|1973 Holweger et al.| 
|1978t |Solanki|[l986| ). Quiet Sun Hinode SOT/SP observations 
containing only weak magnetic fields ha ve been fitted withou t 
the inclusion of a micro turbulence term ( |Socas-Navarro||201 1| ). 
In areas with kG magnetic fields, however, an excess of turbulent 
velocities was found by |Solanki| ( |1986 ). The total rms value for 
turbulent velocities has a range of 1.0 km/s to 3.5 km/s in areas 
containing kG magnetic fields, depending on the spectral line 
and on the line strength. [Zayer et aL] ( |1989| ) found, using lines in 
the infrared, that rms turbulent velocities between 3.0 km/s and 
3.5 km/s were necessary to fit the observed profile shapes. The 
dependence on spectral line strength suggests a height dependent 
turbulent velocity. Therefore, we carried out the inversions with 
a height dependent microturbulence, while forcing the macro¬ 
turbulence to zero. This approach turned out to give satisfactory 
fits to the line profiles. Restricting ourselves to microturbulence 


alone is also in line with the improved spatial resolution and sta¬ 
ble observing conditions of the Hinode satellite. Any remaining 
non-thermal broadening in the spectral line profiles is assumed 
to be caused by unresolved velocities within the resolution ele¬ 
ment. A mixture of up- and downflows with a large correlation 
length along the LOS is less likely at high resolution, than in the 
low resolution data analysed in the earlier investigations. 

Figure ^3 \ displays the microturbulent velocities, ^mic, retrieved 
by the inversion at the three log(T) heights. In all three panels 
areas with an increased are found to coincide with areas 
of strong magnetic fields. The quietest areas are found to have 


the lowest turbulent velocities, supporting the finding of Socas- 
Navarro| ( [2Qll| ) that turbulent velocities are low in the quiet Sun. 


In particular. Fig. [23p reveals that although core pixels display 
enhanced microturbulence, the largest turbulent velocities are 
preferentially situated at the edges of a magnetic feature, form¬ 
ing a narrow halo or ring around MFCs. A closer inspection of 
Fig. [2^ reveals that the same is also true at log(T) = -2.3, al¬ 
though the halo is less well marked at this height. At log(T) = 0 
the halo of large microturbulence around the MFCs is more ex¬ 
tended, although large microturbulent velocities are found every¬ 
where. In this inversion the upper limit for the microturbulence 
was set at 5 km/s. If a higher upper limit was set then microtur¬ 
bulent velocities of up to 9 km/s were retrieved in the halo by the 
inversion at log(T) = 0. 

The distribution of ^mic, found in core pixels, are given by the 
solid histograms in Fig.[^ These distributions can be compared 
to dotted histograms depicting obtained from quiet Sun pix¬ 
els. The differences in between the quiet Sun and MFCs, 
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Fig. 23. a-c: Microturbulent velocities retrieved by the inversion at log(T) = 0, -0.9 and -2.3 from left to right. Note the different 
colour scale for each frame. The black contour encompasses core pixels in all images. The dotted lines (in panel c) display canopy 
pixels. 
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Fig. 24. Histograms of microturbulent velocities in core pixels, solid and the quiet Sun, dotted, at log(T) = 0, a, at log(T) = -0.9, 
h, and at log(T) = -2.3, c. The dashed histogram in panel, c at log(T) = -2.3, represents the microturbulent velocities of canopy 
pixels. 



Microturbulence [km/s] 

Fig. 25. Histograms of microturbulent velocities in pixels imme¬ 
diately surrounding core pixels, same as Fig. [m at log(T) = 0 in 
red and log(T) = -0.9 in green. The dotted histograms show the 
microturbulent velocities in the quiet Sun for the same layers. 


suggested in Fig. are confirmed by Fig. At log(T) = 0 
^mic in the quiet Sun is on average slightly larger, with a mean 


of 3.1 km/s, than in MFCs, which have an average velocity 
of 2.8 km/s. In higher layers the micro- turbulence in MFCs is 
significantly larger than in the quiet Sun. Within the MFCs the 
average ^ynic are 1.9 km/s and 1.3 km/s at log(T) = -0.9 and 
log(T) = -2.3, whilst in the quiet Sun the average ^rnic values are 
0.8 km/s and 0.4 km/s in the two layers, respectively. Note that 
the canopy pixels require a similar microturbulence as the core 
pixels at the same optical depth. Fig. 7^ also reveals that in the 
upper two layers there are many pixels in the quiet Sun which 
require no microturbulent broadening at all. 

The microturbulent velocities found in pixels, which are imme¬ 
diately adjacent to core pixels are displayed in Fig. 25 The pix¬ 
els used in this figure are identical to the ones used in Fig. 
and are, therefore, the pixels featuring the highest LOS down¬ 
flow velocities. Again, only the microturbulence in the lower two 
log{T) layers are shown in Fig. (due to the overlying canopy 
at log(T) = -2.3). The average and median microturbulent ve¬ 
locities are 4.4 km/s and 5.0 km/s at log(T) = 0 and 2.2 km/s and 
2.3 km/s at log(T) = -0.9. These values are higher than those in 
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the quiet Sun and also than those found within the core pixels, 
demonstrating that the highest microturbulence is located at the 
edges of MFCs. Both the average and median values of at 
log(T) = 0 in the surrounding, downflowing pixels are likely to 
be lower limits due to the upper bound of 5 km/s set on the mi¬ 
croturbulence in the inversion. 

The microturbulence retrieved by the inversion could, but does 
not necessarily imply the existence of turbulence or unresolved 
convective processes taking place within the magnetic features. 
It may point to unresolved waves in the MFC (e.g. surface waves 
could account for the higher near the boundaries of the 
MFCs), or it may be due to unresolved horizontal velocity gradi¬ 
ents strongest at the boundaries. Another possibility may be, at 
least in part, a signal of magnetic reconnection between the op¬ 
posite polarity fields found in this study (see Sect. 4.9). Although 
it cannot be completely ruled out that inaccu racies in the damp 
ing constants of Fe I 6301.5 A and 6302.5 A (Anstee & O’Mara 


1995t|Barklem & Q’Mara|1997t|Barklem et al.|1998|) may con 

tribute to the deduced such inaccuracies are unable to ex¬ 
plain the excess in in the magnetic features, since the spec¬ 
tral lines are significantly weakened there, so that damping be¬ 
comes less important. 


5. Discussion 

In the preceding sections we sought to ascertain some of the 
characteristic properties of the magnetic field in solar plage. 
Stokes /, 2, 1] and V profiles of AR 10953 observed by the 
SOT/SP aboard Hinode were analy sed using th e SPINOR 
inversion code ( [Frutiger et al. [2000 ), extended by |van Noort 


( |2012| ), to perform a spatially coupled 2D inversion, which 


Wiehr| (e.g. 1978]); Zayer et al.| 

1990 

); IRiiedi et al. 

(e.g. |1992|( 

1995 

); |Martmez Pi 

let et al.| ( 


(e.g. 1991) . These values are. 


however, somewhat larger th an those found using th e Fe I 
5250.2/5247.1 A line pair (e.g. Stenfio & Harvey[|l985| , which 
is likely related to differences in the heights of formation of the 
two line pairs. Since no magnetic filling factor is introduced the 
obtained B values are lower limits to the true field strengths. 
The magnetic field is found to drop more slowly with height 
when compared to a thin-tube (zeroth order) model. This lack of 
agreement stems from a lack of flux conservation with height, 
by as much as 20%, that primarily affects the upper log(T) node 
in the inversion. If this increase in magnetic flux is compensated 
then the stratification of the strongest magnetic fields with 
strengths higher than 2000 G, found at the centres of MFCs, do 
agree with the thin flux-tube model. When only the lower two 
log(T) layers are considered the model and the inversion also 
agree well. For magnetic fields < 2000 G at log(T) = 0 other 
effects start to play a role. 

The expansion of five isolated MFCs was compared to the 
ideal expansion of two order thin flux tube models. The 
expansion of the selected MFCs and t he models agreed well , 
supporting previous results reported by Pietarila et al. ( |2010 ). 


At the highest inverted layers the field did expand somewhat 
less rapidly than in the model, however, probably due to 
the interactions (merging) with other magnetic features. The 
inversion demonstrates that the majority of MFCs in the strong 
plage investigated here merge with neighbouring MFCs already 
in the middle photosphere. Above the merging height B drops 
more slowl y or not at all and the field also does n ot expand 
very much ( [Pneuman et al.||l986} [Steiner et al.||1986| ). The fact 
that a zeroth order flux tube i s too simple to describe a MFC 
( Yelles Chaouche et al. 2009 ) and expands more rapidly than 
more realistic models that include curvature forces may also 
contribute. 

The LOS velocities obtained from the inversion show that the 
bulk of the magnetized gas in MFCs is essentially at rest and 
shows only weak downflows typically on the order of 200 m/s at 
log(T) = -0.9, which agree s with the results of |SolankI ( |1986| ); 
Martinez Pillet et al. \\991) . Some of the core and canopy pixels 


show upflows or downflows of up to 1 km/s, which were also 
observed by [Langangen et ak ( 2007| ). They are in most cases 


removes the influence of the instrument’s PSF in parallel with 
inverting the data. 

The inversion was able to retrieve and reproduce many of the 
previously deduced characteristic properties of magnetic flux 
concentrations (MFCs) in plage regions. The typical magnetic 
field strength of MFCs was found to be in the kC range and 
took on an average value of 1520 G at log(T) = -0.9, where the 
6302 A line pair is most sensitive to the magnetic field. Similar 
photospheric magnetic field values were found previously by 


found within the MFCs and could be the result of oscillations 
or other transient events occurring within MFCs. The MFCs are 
each surrounded by a ring of strong downflows, which can be 
readily observed at the log(T) = 0 and -0.9 nodes in Figs. 

These downflow rings were also seen to shift outwards in 
higher log(T) layers as the MFC expands. A downflow ring has 
an average velocity of 2.4 km/s at log(T) = 0, but parts of the 
ring achieve supersonic velocities up to 10 km/s, corresponding 
to a Mach number of 1.25, in the same log(T) layer. On average 
2.5% of a MFC’s ring contains downflows reaching supersonic 
velocities. Some of these downflows appear to overlap with the 
magnetic field, which could be the signature of entrainment, 
or may be a result of insufficient spatial resolution to cleanly 
resolve the magnetic boundary of MFCs. It is unlikely that 
they indicate convective collapse, since the Stokes V profiles 
of these p ixels do not have th e characteristic third lobe re¬ 
ported by [Nagata et al.| (|2008|). Downflows at the edges of 


magnetic features were previously found by |Rouppe van der 

|Voort et al.|(|20051) and Langangen et al.|(|2007 

). They are also 
le modelling of 

compatible with the downflows inferred from t 

Stokes y asymmetries (|Grossmann-Doerth et a] 

l.|1988HSolanki 


log(T) = 0, gradually weaken higher up in the atmosphere 
and have completely disappeared at log(T) = -2.3, suggesting 
that they are associated with the granular convection pattern. 
However, downflows exceeding 6 km/s as well as supersonic 
velocities are only found at the boundaries of MFCs and never 
in the quiet Sun, indicating that the granular downflows are 
strengthened by the presence of the MFC. We may be seeing 
a similar effect as at the edges of light-bridge granules, which 
display extremely strong downflows, probably because of a 
combination of radiative losses into the magnetic feature, which 
leads to faster flows, and because the Wilson depression allows 
us to see deeper layers with faster flows ( |Lagg et al.|2014| ). The 
downflows around MFCs also appear qua litatively similar to 
the downflows observed around pores by [Hirzberg^ ([2Q03|); 
[S^ankarasubramanian & Rimmele] ( |2003] ) and |Cho et al.| ( |201^ 


and, therefore, might show a high temporal variation as well. 
Several MFCs also possess a plume-like downflow features 
characterised by localised strong downflows traceable through 
all log(T) layers. These plumes become larger in size and 
smaller in magnitude with height. They also trace the expansion 
of their MFC by shifting outwards with height. The origin of 
these plumes might be tied to a particular active granule in the 
vicinity, or may be caused by plasma pouring down from the 


17 











































































































D. Buehler, et al.: Properties of solar plage 


chromosphere. A reconnection event in the chromosphere might 
also give rise to such a structure, but the rapid increase of LOS 
velocities with depth appears to be incompatible with such a 
scenario. 

Located within the downflow rings surrounding the MFCs, we 
found small magnetic patches bearing the opposite polarity to 
the main MFC to which they adjoin. These opposite polarity 
patches are only visible in the lower two log(T) nodes, as 
shown in Fig.[^ supporting the notion that they are intimately 
connected to the downflows in which they are immersed. MHD 
simulations carried out by Steiner et al. ( |1998| ) and |Vogler et al.| 
( |2005| ) predict such patches in th e vicinity of strong magnetic 


field concentrations, (see Fig. 5 in Vogler et~ST]|2005|). Indirect 


could be identified for all MFCs as well as for the pores and 
sunspot. The canopies were found to harbour weaker, more 
horizontal magnetic fields, with inclinations as high as F = 80°. 
The typical LOS velocities in the canopies are identical to the 
LOS velocities found at the core of a MFC. Many magnetic 
canopies were found to lie above essentially field free regions. 
The canopies found here agre e with the model proposed by 
Grossmann-Doerth et al.| ( T958| ) for the production of the Stokes 


V area asymmetry and the observational results pre sented by 
Rezaei et n^(|2007|);|Narayan & Scharm^ ( |2010| ) and [Martini 


Gonzalez et al.| ( 2012| ). All magnetic features were observed to 


evidence was also found in a network patch by |Zayer et al.| 
( [1989| ), although the exact spatial proximity to the MFC could 
not be established in that publication due to the low spatial 
resolution. However, the reversal of the Stokes V amplitude in 
the wings of the 630 nm line pair is only seen in the deconvolved 
profiles produced by the inversion, where the Stokes V spectrum 
displays three lobes. An additional complication is the masking 
of these magnetic fields by the canopy of the main MFC. 
The canopy has a typical field strength of 300 G or higher at 
locations overlying such weak opposite polarity fields. The 
canopy thus produces a strong signal in Stokes Q, U and V, 
whereas the field strengths of the opposite polarity patches at 
log(T) = 0 are typically well below canopy values. This also 
prevents the selection of these patches via a typical amplitude 
threshold in the Stokes profile, so that we cannot completely 
rule out that even those opposite polarity patches with 5^100 G 
are an artifact of the inversion. The opposite polarities typically 
carry less than 1% of the fiux present in their parent MFCs. 
Further observations performed with a higher spatial resolution 
are necessary to ascertain the existence of these small opposite 
polarity patches. Such observations would be particularly useful 
in the 1.56 pm lines, due to their large Zeeman se nsitivity and 
low fo rmation height. The 1.56 pm data analysed by |Zayer et al.| 
( |1989| ), although of low spatial resolution, provide some support 
for our results. The close coexistence of such opposite polarities 
might give rise to current sheets and could lead to reconnection 
events. 

We have introduced a novel method for the resolution of the 
180° azimuth ambiguity applicable to largely unipolar regions. 
It makes use of the basic thin fiux-tube structure of MFCs and 
assumes that the divergence of the magnetic field tends towards 
zero. With this method we find that the core pixels of MFCs 
have typical inclinations, relative to local solar coordinates, 
between 10° and 15° in all three log(T) nodes of the inversion, 
which agrees well with earlier inclination results found by 
Topka et al. (|1992| ); |Bernasconi et al.| ( |1995 ) and |Martmez| 
Fillet et aL] ~( 1997| ) The distribution of the azimuths shows a 
preference for the eastern direction, which can be attributed to a 
LOS effect. MFCs located closest to the disc centre in the field 
of view showed the most homog eneous azimuth distribution, 
which supports the conclusion of Martinez Fillet et al. ( \991) 
that in general MFCs have no preferred orientation. MFCs close 
to the sunspot did, however, display azimuth distributions that 
were either predominantly directed towards or away from the 
spot depending on their polarity. The canopies of these MFCs 
were also irregular and elongated. This demonstrates that a 
nearby sunspot has a direct impact on the properties of MFCs 
even in the middle photosphere. 

The inversion allowed us to clearly differentiate between 
magnetic fields which form magnetic canopies from those 
which form the core or root of a MFC. A separate canopy 


expand with height and many isolated MFCs merged to form 
comparatively large expanses of magnetic field at log(T) = -2.3. 
We therefore expect that at least some of the magnetic features 
display a similarity to the wine-glass model described by|Bunte 
[etaL] ( [T99^ . 

The average temperature stratification within MFCs most 


closely followed the empirical plage fiux-tube model of lSolanki 
& Brigljevic ( 1992| ). Temperatures were hotter by 300 K within 


the MFCs when compared to the quiet Sun in the upper two 
log(T) layers. However, a closer inspection of MFCs revealed 
that the temperature is not homogeneous in any log(T) across 
a MFC. The highest temperatures in a MFC are typically 
positioned at its edges and concentrated into isolated points at 
log(T) = 0. Often supersonic downflows were located in the 
vicinity of these temperature enhancements. This raises the 
possibility that the two might be connected, since the supersonic 
flows likely create a shockfront below the log(T) = 0 layer. 
However, there were also cases of temperature enhancements 
without a nearby supersonic velocity. It is likely that a clear 
picture of the effect of these supersonic downflows upon an 
MFC can only be determined together with an analysis of a time 
series of data. 

The inversion procedure included a depth-dependent mi- 
croturbulent velocity, with the macroturbulence set to zero. 
Farticularly high values of the microturbulence, partly in 
excess of 5 km/s at log(T) = 0, were found at the edges of 
the magnetic features. The typical microturbulent velocities in 
MFCs obtained by the inversion has an average value of 2.8 
km/s at log(T) = 0 and becomes weaker with height. Such 
values are on a par with the br oadenings presented by |Solanki| 
( |1986| ) and Zayer et al. ( |1989| ). This suggests the presence of 
strong unresolved velocities around the MFCs. Also, although 
its amplitude decreases rapidly with height, ^mic is at least 
a factor of 2 larger in the MFC than in the quiet Sun in the 
upper two log(T) layers. These large microturbulence values 
may be telling us that the strong downflows at the edges of 
the MFCs are associated with vigorous turbulent motions, or it 
may be a signature of waves travelling along the MFC. Finally, 
the microturbulence returned by the inversion code may be 
signalling strong horizontal velocity gradients, e.g. across the 
boundary of MFCs. A possible geometry for such gradients is a 
strong downflow outside and a weak one inside. 


6. Conclusion 

In this investigation we have been able to confirm many previ¬ 
ously obtained properties associated with magnetic fields, which 
form solar plage. At the same time, we have also uncovered new 
features associated with kG magnetic fields. KiloGauss fields 
in plage are strong, vertical (10° - 15°) magnetic fields on the 
order of 1.5 kG, which expand with height similar to a thin 
fiux-tube forming extensive canopies as low as the upper photo- 
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sphere. These fields are further characterised by being on aver¬ 
age 300 K hotter than their surroundings in the middle and upper 
photosphere, in agreement with empirical fiux-tube models, and 
contain weak, on average 200 m/s, plasma fiows within them. 
This investigation has also discovered several new properties of 
these fields. Each magnetic fiux concentration is surrounded by 
a ring containing strong downflows, on average 2.4 km/s in the 
lower photosphere. A typical ring shifts outward with height as 
the field expands and parts of it can attain supersonic veloci¬ 
ties in the lower photosphere. Co-spatial with these rings we 
found enhanced microturbulent velocities decreasing with height 
as well as magnetic patches situated beneath the canopy with a 
polarity opposite to the main plage-forming fields. The plasma 
temperature of the MFCs was not uniform across their cross 
sections, displaying instead temperature enhancements at their 
edges. Magnetic field concentrations located close to larger fea¬ 
tures such as pores or a sunspot displayed an asymmetric canopy 
and more inclined magnetic fields. 
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